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Abstract 

We characterize the multipartite entanglement in a quantum system by the quantity 
which vanishes if only the quantum system may be decomposed into two weakly entan- 
gled subsystems, unlike measures of multipartite entanglement introduced before. We 
refer to this quantity as the minimal entanglement of bipartite decompositions (MEBD). 
Big MEBD means that the system may not be decomposed into two weakly entangled 
subsystems. MEBD allows one to define, for instance, whether the given quantum system 
may be a candidate for a quantum register, where the above decomposition is undesir- 



C^ I able. A method of lower estimation of MEBD is represented. Examples of big MEBD in 

spin- 1/2 chains governed by the Hdz Hamiltonian in the strong external magnetic field 
are given. 



^. 1 Introduction 

-(— > ■ 

Ch ■ The problem of strong entanglement in spin systems is very important in view of development 

S . of quantum computation. It is acknowledged that the quantum correlations are needed in order 

O^! to organize the coherent manipulations by different bits in quantum circuit. This is a basic 

I resource of quantum computation providing advantages of quantum circuits in comparison 

^ ' with their classical counterparts [1]. In particular, quantum correlations are responsible for 

^^ • the speedup of quantum algorithms. It is hard to create large strongly entangled systems. 

cn . On the contrary, a few-qubit quantum registers have been constructed using different physical 

(^ ! basis. For instance, there are registers using systems of coupled trapped ions [2], registers 

ly-T I using superconducting charge qubits [3] and neutral atom quantum registers [Ij. These small 

O ] registers can be organized in large scale quantum circuits [5], El [7] . Review of different measures 

2 ' of both bipartite and multipartite entanglements is given, for instance, in [8],|9]. Especially one 

has to mention the Wootters criterion which is well elaborated for calculation of entanglement 

among two nodes in a spin- 1/2 system [lOl [11]. The entanglement between two subsystems 

rS I consisting of more then one node may be effectively described by the positive partial transpose 

c3 ; (PPT) criterion [121 IE] introducing so-called double negativity as a measure of entanglement. 

This measure is applicable to both pure and mixed states. An important concept is so-called 

entanglement of formation [H] describing the entanglement between two subsystems in a mixed 

state. 

This paper is devoted to the problem of multipartite entanglement. Whereas the bipartite 
entanglement is studied in many details, the multipartite entanglement is more cumbersome 
[T^ [TB] and it has not been completely understood yet. Of course, considering the problem 
of strong multipartite entanglement in a quantum system, it would be reasonable to construct 
such strongly entangled system, where any two nodes are strongly entangled simultaneously. 
However, this requirement seems to be too tough. Therefore we assume that the strong node- 
to-node entanglement between any two nodes is not necessary in order to entangle all nodes in 



the quantum system. Instead of this, we state that the system becomes strongly entangled at 
some time moment to if any its subsystem is strongly entangled with the rest of the system at 
this time moment. 

We suggest a method to pick out such quantum systems which may not be decomposed into 
two weakly entangled subsystems. It is shown that such systems may be simply realized using 
spin chains. Entanglement considered here is some kind of "collective" entanglement because 
not all node-to-node entanglements in such systems are big. We assume that such "collective" 
entanglement may provide advantages of quantum devices in comparison with they classical 
counterparts. 

Following the above two paragraphs, we introduce so-called minimal entanglement of bipar- 
tite decompositions (MEBD) of the quantum system as a witness of strong entanglement, or 
as the measure of the above mentioned "collective" entanglement. The basic feature of MEBD 
is that its measure vanishes if only there is at least one possibility to decompose the given 
quantum system into two weakly entangled subsystems. Thus, MEBD helps us to test, for 
instance, whether the given quantum system may be a candidate for a large volume quantum 
register, where such decomposition is undesirable since one has to be able to address to all 
qbits of register at ones with a single operation. Namely this property of quantum register is 
responsible for the exponential speedup of quantum algorithms in comparison with their classi- 
cal counterparts. In other words, if MEBD vanishes, then the system must be considered as two 
(or even more) weakly entangled subsystems rather then a single entangled quantum system. 
Emphasize that big MEBD does not require significant entanglement between any two particu- 
lar nodes simultaneously. Instead of this, a big MEBD means that the whole quantum system 
is entangled and there is no decomposition of this system into two un-entangled subsystems. 

We consider the dynamics of MEBD in homogeneous spin-1/2 chains governed by the H^z 
Hamiltonian in the strong external magnetic field. It will be demonstrated that big MEBD 
is achievable during the relatively short time interval (comparable with the end-to-end state 
transfer time interval r along the properly adjusted inhomogeneous spin-1/2 chain, r ~ vr 
[171 [18]) and this time interval is slightly increasing with the length of the spin chain. It is also 
remarkable that big MEBD is achievable in relatively simple models such as the homogeneous 
spin-1/2 chains with special initial conditions. These special initial conditions are represented 
by the certain number of the initially excited spins, i.e. by the number of spins which are 
directed opposite to the external magnetic field initially. Due to these initial conditions we are 
not forced to use either the inhomogeneous chains [191 120] or inhomogeneous magnetic field [21] 
to obtain big MEBD, unlike the node-to-node entanglement in the spin chains with the single 
excited node which has been considered, for instance, in [T^ I2U1 1^ . 

The definition of MEBD is introduced in SecUl An algorithm allowing one to obtain a lower 
estimations of MEBD is given in Sec J2.1[ This lower estimation becomes important for large 
systems, where the calculation of MEBD following its definition is very complicated. Examples 
of big MEBD in homogeneous spin-1/2 chains governed by the Hdz Hamiltonian in the strong 
external magnetic field are represented in SeclHl We restrict ourselves to the short spin chains 
because even multipartite entanglement in the short chains is not a simple phenomenon. For 
instance, it is impossible to organize the strong entanglement between any two nodes in the 
three node spin-1/2 chain simultaneously, while big MEBD is achievable in three- node and 
larger spin chains. Basic results and conclusions are summarized in SecJH 



2 The minimal entanglement of bipartite decompositions 
of a quantum system 

Hereafter we will characterize the entanglement between two subsystems a and h of the iV-node 
quantum system Sjq = a U 6 (a fl 6 = 0) in terms of the double negativity Na,h = ■^h,a (PPT 
criterion [111 [13]), where Nafi is the absolute value of the double sum of the negative eigenvalues 
of the matrix p-^" (or p^*) which is the transposition of the density matrix p (associated with 
the system Sn) with respect to the subsystem a (or h). More general, we will also use the 
bipartite double negativity Ma^^a- in the L-partite system S'tv, 

S'at = ai U 02 U ■ ■ ■ U a^, a^ fl aj =0, if z 7^ j. (1) 

In this case A/'a-^a = -^a^ai is the absolute value of the double sum of the negative eigenvalues of 

the matrix p^^^ (or Pj/' )• Here pij is the density matrix p reduced with respect to all subsystems 
Cfc except tti and af pij = Ti {ak,kj^i,j}P ^^^ superscript T^. (or T^J means the transposition 
with respect to the subsystem aj (or Oi). The choice of PPT criterion provides applicability of 
the represented algorithm to both pure and mixed states. 

The results of this paper are based on the following definition. 

Definition: Let us consider the set of bipartite entanglements of all possible bipartite 
decompositions of the given quantum system. The minimal entanglement out of this set will 
be referred to as the minimal entanglement of bipartite decompositions (MEBD) of the given 
quantum system. 

In accordance with this definition, the big value of MEBD means that the quantum system 
may not be decomposed into two weakly entangled subsystems. 

To clarify the importance of MEBD we refer to the quantum register. Remember, that the 
inherent property of quantum register is the simultaneous entanglement among all its nodes. 
This entanglement allows one to reach the exponential speedup of the quantum calculations in 
comparison with the classical ones. Thus, if the A^-node quantum system may be separated 
into two weakly entangled subsystems of A''! and A'^2 nodes (A''i + A'^2 = N), then the effective 
volume of the register constructed on the basis of this quatum system is either Ni or A'^2 rather 
then A^. By definition of MEBD we see that such separation is impossible if only MEBD is big. 
Thus, MEBD helps us, for instance, to test whether all nodes of the quantum system may be 
effectively used as different bits in a quantum register. 

In accordance with the above definition of MEBD, we introduce a measure of MEBD which 
vanishes if only there is at least one decomposition of the quantum system into two weakly 
entangled subsystems (note that the measure of multipartite entanglement introduced in refs. 
[15], [16] does not have this property). Namely, let the system Sn have M different decomposi- 
tions into two subsystems: 

SM = A^'>UAf\ Af^Af =0, t = l,...,M, (2) 

^_iE^lVC'N + lCr = 2^-^-1, iV = 2,4,6,... 
I EE7^^/'C7^ = 2^-^-1, iV = 3,5,7,... ' 

j! 
where C!- are binomial coefficients: C] = ——. -7. Then we introduce the measure of MEBD 



E{Sn) by the following formula: 



^!(J 



E{Sn) = . min AT.a) .(2). (3) 



i=l,...,M ^-- 



This formula is consistent with the definition of MEBD, because the measure is zero if only 
there are at least two un-entangled subsystems AlJ and AiJ: -Afji) A2) = 0. We consider that 

the value of MEBD is big if its measure approaches unity, and the value of MEBD is small if 
its measure tends to zero. 

Emphasize that the big MEBD does not require the strong entanglement between any two 
nodes of the quantum system. In other words, MEBD may be big even if some particular nodes 
are weakly entangled between each other. This statement may be simply justified using the 
property of the hierarchy of double negativities [13], which can be formulated as follows. Let 
the quantum system be separated into the set of subsystems ([H). Then 

J^ai,{a2,...,L} > ■^fai,{a2,...,L-l} > " " " > ■^ai,a2- (4) 

Let us apply this property to two simple examples of spin-1/2 systems with big MEBD 
where only some two-node double negativities are big. Hereafter we will use notation S^ = 
{si, . . . , sn} for the A^-node spin system, where Sj means the zth node. 

First of all, we turn to the three- node spin system 5*3 = {si, S2, S3}, which may be decom- 
posed as follows: 

83 = SiU {S2, S3} = S2 U {Si, S3} = S3 U {Si, S2}. (5) 

In order to obtain big MEBD defined by Eq.Q it is enough to have two big double negativities: 
Afsi,s2 ^-iid J\fs3,{si,s2}- III ia.ct, using the property of the hierarchy of double negativities (J4]) we 
obtain that A/^i, {52,^3} > ^fs^,s2, Ma, {51,^3} > ■f^s2,si = MsuS2- Thus, all J^s„{sj,sk} are big (all 
i, j and k are different) and the system may not be decomposed into two weakly entangled 
subsystems, i.e. MEBD of the system 5*3 is big. 

The second example regards the four node system 6*4 = {si, S2, S3, S4}, which may be de- 
composed into two subsystems as follows: 

5'4 = Si U {S2, S3, S4} = S2 U {Si, S3, S4} = S3 U {Si, S2, S4} = (6) 

S4 U {Si, S2, S3} = {Si, S2} U {S3, S4} = {Si, S3} U {S2, S4} = 
{S1,S4}U{S2,S3}. 

In order to achieve big MEBD, it is enough to have three big double negativities: Msi,s2y ■Mi3,s4 
and M^si,s2}, {33,34}- Then, using the property of the hierarchy of double negativities (j4]), one can 
show that the double negativities for all decompositions ([6]) are big. In fact: 

-'^Si,{s2, •53,54} — -'^Sl,S2) ■'^S2,{si,S3,S4} — ■'^32,31 ^-'^Sl,S2! \' ) 

■'^ S3, {31,32,34} ^ -/^ 33,341 ■'^ 34, {31,32,33} ^■'^34,33 = -/^ 33,341 
■'^{31,33}, {32, 34} ^■'^31,321 ■'^{31, 34}, {32,33} ^■'^Sl,32- 

In result, since J\fsi,s2i -^33,34 and J\f{si,32},{33,34} are big by our assumption, then the above 
system S4 may not be decomposed into two weakly entangled subsystems, i.e. MEBD of the 
system 5*4 is big. 

2.1 Lower estimation of MEBD 

Note that E{Sn) given by formula ([3]) is cumbersome for calculations if A^ is big, because 
M increases very quickly with the increase in A^. Thus, it is reasonable to introduce another 



quantity giving a lower estimation of MEBD in the system Sn, which is based on the following 
remark. 

Let us fix some decomposition of S^, say, S^ = A\ U A\ . Let E{A^ ), i = 1,2, 
be measures of MEBD of subsystems A^ (introduced in the previous decomposition), see 
Eq.(l3]) with replacement Sn by A^ . Then we state that if the quantity E^ (Sn), El (Sn) = 
mm{E{A'{^), E{Af'^),Af,(i) ,(2)), is big, then MEBD of Sn is big, i.e. E{Sn) is big as well. The 

opposite is not true in general, i.e. if E[ (Sn) is small, then MEBD of the system may be 

big. This happens if there is another decomposition Sn = A2 U A2 such that £"2 {Sn) = 

lainiEiAo ),E(A2 ),A/'.(i) .(2)) is big. In other words we may estimate MEBD of the system 

^2 '^2 
of A^ nodes as follows. 

Let us fix the jth decomposition of Sn into two subsystems, Sn = Aj U Aj , I < j < M, 

and consider the following quantities: 



^« 



(Sk) = min {E{Af>), E{Af')M^ai^m) , (8) 



where E{A^ ) are measures of MEBD of subsystems A^ . The manifold of all possible E- (Sn) 
will be referred to as S'^^^Sn)'- £^^\Sn) = {Ej (Sn), J = 1, • • • , M}. It is clear that the value 
E^^\Sn), calculated by the formula 

E^'\Sn)= max e\'\Sn), (9) 

j=l,...,M 

satisfies the inequality E^^^Sn) < E{Sn) due to the hierarchy of double negativities (jl]). Thus 
E^^\Sn) may be referred to as the lower estimation of MEBD of the system Sn- 

We may introduce another lower estimation as follows. Let us decompose each subsystem 
A- and A, into two smaller ones: 

4' = 4!' U A^', A^> n A^ = 0. 1 < * < M<" (10) 

Af = ^™u4f, ^™n^<f =e, l<n<Aff. (11) 

Here MJ , Z = 1,2, are the numbers of all possible decompositions fITOl) and flTTl) respectively 

with fixed j. Let E^, {Aj ) be lower estimation of measure of MEBD for the subsystem A^ , see 
Eq.® for the definition of this estimation. Then another lower estimation may be introduced 
by the following formula: 

E^liSN) = min (i?«(A«), E(i)(Af ), Ar^(i)_^(2,) . (12) 

The manifold of all possible E^J^ISn) will be referred to as S^^\Sn): S^^\Sn) = {E^Z{Sn), 1 < 
j <M, l<k< Mf\ 1 < n < Mf^}. We define the quantity E^'^\Sn), 

E^''\SN)=max Ef^^iSN). (13) 

j,k,n ■' 

Again one has E^'^\Sn) < E^^\Sn) < E{Sn) due to the hierarchy of double negativities (jl]). 
Thus, the quantity E^'^^Sn) is another lower estimation of MEBD of the system Sn- This pro- 
cess may be continued involving further decompositions of the subsystems Ajl"^' into the smaller 



ones to find the manifolds £^p\Sn) and appropriate minima E^p\Sn), P = 3, 4, . . . , Jmax, wliere 
^(Jma3:)(^^^~) corresponds to tlie case wlien all subsystems consist of two nodes. Again we obtain 
tliat 

^(J™.)(^^) < E^-^— -i)(^jv) < ■ ■ ■ < E'^^^Sn) < E{Sn) (14) 

due to the hierarchy of double negativities (jll). 

All in all, if there is a big element eo G £^''\Sn), then we may write cq < E^'^^Sn) < E{Sn), 
i.e. Co may serve as the lower estimation of MEBD of our system. This low estimation will be 
found in examples of SecJS] 

It is worthwhile to clarify, whether the values E^^^'^Sn) can be used for the reasonable 
approximation of E{Sn) or their are just lower bounds of MEBD. First of all, it is necessary to 
remark that E^'''>{Sn) is significantly less then E{Sj\i), which is a consequence of the additivity 
property of the entanglement [23]. On the other hand, this additivity property allows us to 
improve estimation as follows. For simplicity, we consider the case k = 1 and symmetrical 
decomposition Sj^ = A[' U A[' . Of course we may not apply the additivity property to our 

case directly because, for instance, p ^ p^^ ® p^^ . Nevertheless we use it for the rough 
estimations. Assume that M a\) ai) > E{A\ ), which is reasonable for the big MEBD when 

E{Sn) ~ 1. Then E^ (Sn) = E{A[ ). In accordance with the additivity property we may 
write 

E{Sn) ~ 2E(4')) => E(4')) ~ ^E{Sn). (15) 

Consequently we may propose that E^ (S^) (and E'^^^Sn)) is big if Ei (Sn) ~ 1/2. In this 
case E{Sn) ~ 1. Thus E^ (Sn) may be used to approximate E{Sn). 

3 The dynamics of MEBD in the homogeneous spin- 1/2 
chains governed by the Hdz Hamiltonian. 

The dynamics of entanglement significantly depends on the physical scenario, in particular, 
on the initial state of the quantum system. For instance, the dissipative dinamics of multi- 
partite entanglement is considered in set of papers, see [22l [23]. In this section we consider 
the quasiperiodic dynamics of MEBD caused by the special initial condition in the spin- 1/2 
homogeneous chain governed by the H^z Hamiltonian 

^ 2fc 

in the strong external magnetic field. Here 7 is the gyromagnetic ratio, rij are distances 
between the zth and the j'th nodes, h is the Planck constant. Hamiltonian Hdz is the secular 
part of the Hamiltonian describing the dipole-dipole interaction in the strong external magnetic 
field [25]. This Hamiltonian is known to commute with Iz (the ^-projection of the total spin). 
The density matrix, obtained as the solution to the Liouville evolution equation ipt = [Hdz, p] 
{h= 1), reads 

pit) = e-'^-'^' p^e^-^\ (17) 



where po is the initial density matrix. In the homogeneous chain -Di,i+i = D, i = 1, . . . ,N — 1. 
Hereafter we use the dimensionless time r = Dt. It is pointed above, that the dynamics of 
entanglement (and, consequently, the dynamics of MEBD) strongly depends on the initial state 
of the spin chain, in particular, on the number of initially exited spins (i.e. the spins which are 
directed opposite to the external magnetic fields). We use notation |\l/o(S'Ar)) for the initial 
wave function, so that po{Siy) = \'^q{Sn)){'^o{Sn)\- We have found that in order to achieve 
big MEBD at some time moment in the homogeneous spin chains with A^ = 3,4,6,8 one has 
to take the initial condition |\l/o(S'Ar)) with one, two, three and four exited spins respectively. 
We consider that MEBD is big at some time moment tn if the introduced measure E{Sn) 
(see Eq.(l3])) approaches unity at this time moment. Emphasize, that the problem of the strong 
simultaneous entanglement between any two nodes is not resolved even for the spin chains of 
several nodes. For this reason the observation of big MEBD in the short chains is valuable. It is 
found that a big MEBD may be achieved even though the entanglement between two particular 
nodes is weak. In all examples given below, tn, N = 3,4,6,8, means the time moment when 
big MEBD is achieved for the first time in the homogeneous A^-node spin chain. Hereafter we 
do not write r as argument of functions for the sake of simplicity. We use the Dirac notations 
\ni . . . un) for the wave functions, where n^ = if the ith spin is directed along the external 
magnetic field, or rij = 1 if the ith spin is directed opposite to the external magnetic field, i.e. 
the zth spin is excited. 

Example 1: A^ = 3. As a simple example we consider the chain of three spins, 53, with the 
initial state |\l/o(5'3)) = |010). All possible decompositions of 83 into two subsystems are given 
in Eq.([5]). The evolution of ^(6*3) = mm{Afs-,^{s2,s3}^-^s2,{si,s:i},J^s:i,{si,s2}) (see Eq.([3])) is shown 
in FiglU^a). The first maximum of this function is achieved at t^ = 1.505, £{3^) = 0.943, i.e. 
MEBD is big at T3. 

Example 2: A^ = 4. Consider the chain of four spins, 5*4, with the initial state |\l/o(5'4)) = 
1 1001). In this case all decompositions of S4 into two subsystems are given in Eq.([n]). Note, 
that the double negativities Msj,{s^,,Si ,Si } (all j,k, 13,14 are different) are less then the double 
negativities associated with other decompositions in set ([6]), except the short interval in the 
vicinity of r = 0. Using definition ([3]) we calculate the evolution of £{84) which is shown 
in Fig{T|^6). The first maximum of this function is achieved at t^ = 1.819, £{84) = 1.000, 
i.e. MEBD is big at T4 so that the system may not be decomposed into two weakly entangled 
subsystem at the time moment T4. 

Let us calculate the lower estimation of MEBD using the technique developed in Sec J2.1[ 
Let A^^^ = {si,S2}, A^^ = {53,54} and calculate ^iH^4) = mm(E(A^^^), E(Af^),J\f,(i) ,(2)) = 

min(A/'si,s25-^s3,s4)-^{si,s2}, {53,54})- The evolution of this function is shown in Fig{T]^6)), dotted 
line. 

Example 3: A^ = 6. Consider the chain of six spins 5*6 with the initial state |\l/o(5'6)) = 
1 100110). The system may be decomposed into two subsystems as follows: 

5*6 = SiU {Sj,Sk,Sl,Sn,Sni} = {Si,Sj}U {Sk,Sl,Sn,Sm} = (18) 

\Si, Sj, Sfcj U j^S;, Sn, Syjij, 

all i,j, k, I, n, m are different, z, j, k,l,n,m = 1, . . . , 6. 

Note, that the double negativities ■f^Si,{s-,Sk,si,s„,sm} i.^^^ i,j,k,l,n,m are different) are less then 
the double negativities associated with other decompositions in set flTSl) . except the short in- 
terval in the vicinity of r = 0. Using definition ([3]) we calculate the evolution of E{Sq) which is 
shown in Fig{T]^c). The first maximum of this function (i.e. a big value of MEBD) is achieved 

7 




Figure 1: The evolution of the entanglement E(Sn) (solid lines), of the lower estimation of the 
entanglement El (Sn) (dotted lines) and of E{Sn) = mini{Afg_,^^g^N) (dashed lines) in chains 
Sn with the special initial conditions |^o(57v)), N = 3,4,6,8: (a) N = 3, |^o(53)) = |010); (6) 
iV = 4, 1^0(54)) = |1001); (c) iV = 6, 1^0(56)) = |100110); (d) N = 8, I'^oiSs)) = |10011001) 



at re = 2.110, EiSe) = 0.992. 

Let us find the lower estimation of MEBD. Let A[' = {si, S2, S3, S4}, A\' = {55, Sg} and 
calculate ^K^e) = min(E(AS^^), E(Af^),7V,(i) ,(2)). Here E{A'i^) is given by eq.® with de- 

compositions of A\ given by eq.(|6]), while E{A[ ') = A/s^e- The evolution of EI{Sq) is shown 
in FigJU^c), dotted line. 

Example 4: N = 8. Consider the spin chain 5*8 with the initial state l^ol'S's)) = 
1 10011001). The system may be decomposed into two subsystems as follows: 

\Si, Sj, Sis J U \Sly Sfiy Sfri) "Spi ^qi \^i) ^ji ^ki ^l J U I'^n; S-m, Sp, Sqj, 

all i,j, k, l,n,m,p, q are different, i,j, k,l,n,m,p,q = 1, . . . , 8. 

Note, that the double negativities ■f^si,{s-,sk,si,s„,s„^,sp,sq} (all i,j,k,l,n,m,p,q are different) are 
less then the double negativities associated with other decompositions in set (TT9|) . except the 
short interval in the vicinity of r = 0. Using definition ([3]) we calculate the evolution of E{Ss) 
which is shown in Figll]((i). The first maximum of this function is achieved at Tg = 2.193, 
E{Ss) = 0.988, i.e. MEBD is big at rg. 

We also calculate El(Ss) (see Eq.®) with A\ = {51,52,53,54} and A[ = {55,56,57,58}. 
The evolution of EUSs) = mm(E(A^^^),E(Af'^),Af,w ,(2)) is shown in FigJITrf), dotted line. 

Here E(A\ ) and E{A\ ) are minima of all double negativities associated with all possible 
bipartite decompositions of A{ and Ai , where the decompositions of Ai coincide with those 
given in Eq.Q, while the decompositions of Ai are given in Eq.Q with replacements 5j — )■ 
5i+4, i = 1,2,3,4. 

Evolutions of MEBD is shown in FigJU (solid lines) demonstrates that the graphs of MEBD 
for systems with N > 3 oscillate around some value in the vicinity of unity. The amplitude of 
these oscillations decreases with the increase in A^ so that MEBD tends to straight line E = 1. 
Regarding the lower estimations (dotted lines), one must note that they do not approach unit. 
They also oscillate and the amplitudes of these oscillations decrease with the increase in A^. 
The low estimation do never approach zero for N = 8 and r > 0. It may be considered as a 
function oscillating around some value Ei ^ 0.4. Comparison of solids and dotted lines shows 
that the maxima of E(Sn) and E^ (Sn) well correlate with each other, although the maxima 

of E(Siy) are less localized then those of Ei {Sjy). The fact that the low estimation is less then 
MEBD is a consequence of the fact that MEBD is a measure of the " collective" entanglement 
which is carried by the dencity matrix p. At the same time, calculating the low estimation, 
we use the reduced dencity matrices which do not contain the complete information about the 
"collective" entanglement. This is in agreement with the estimation given in eq. (fT5!) . 

It is also important to note, that the evolution of MEBD shown in FigH] almost coinsides 
with the evolution of the quantity E{Sn) = m.iia{Mg_ rest'^) (dashed lines), where restf means 

the whole system S^ except the node 5j (solid and dashed lines coinside almost everywhere). 
If so then one can use E{Si\i) in order to estimate MEBD with high precision. However, this 
statement remains unproved for an arbitrary quantum system. 

The first maxima of the functions E{S]\f) together with the appropriate time moments are 
collected in Table 1. It is important to realize how long are the found time intervals tat. 
These intervals must be short enough in order to provide the coherent manipulations by qbits 
in a quantum circuit. For this purpose we note that the dynamics of MEBD is quasiperiodic. 



similar to the dynamics of the quantum state transfer probabihty along the spin chain. Since 
both MEBD and quantum state transfer deal with the same physical object (spin chain) it is 
reasonable to compare the time intervals rjy (obtained in our paper) with the time intervals 
required for the end-to-end quantum state transfer. This quantum process is well studied 
[T7t \TE\ dSl [201 EI] and the minimal possible time interval is established [T7] : it equals to vr in 
dimensionless variables. In our case, t^ < vr, which indicates that the big values of MEBD are 
achievable during the reasonable time intervals. 



A^ 


3 


4 


6 


8 


E{Sn) 


0.943 
1.505 


1.000 
1.819 


0.992 
2.110 


0.988 
2.193 



Table 1: The maximal values E{Sn), N = 3,4, 6, 8, and the appropriate time moments tn 



4 Conclusions. 

We introduce MEBD as a witness of the strong entanglement in a quantum system. Although 
the strong entanglement between any two nodes of a quantum system is hardly achievable 
simultaneously, the big MEBD is quite realizable. The basic feature of MEBD is that its 
measure vanishes if only there is at least one decomposition of the quantum system into two 
weakly entangled subsystems. This feature does not appear in other measures of entanglement. 
Having this property, MEBD may be used, for instance, to test whether A^-node spin system 
may be a candidate for the A^-qbit quantum register or this quantum system must be separated 
into two registers of smaller size. 

We suggest a method of lower estimation of MEBD which is useful for the large systems. 

The dynamics of double negativities in the homogeneous spin-1/2 chains governed by the 
Hdz Hamiltonian with the appropriate initial conditions demonstrates us a possibility to obtain 
big MEBD of a quantum system during the relatively short time intervals tn < vr, which is 
comparable with the shortest time intervals obtained for the quantum state transfer along the 
spin chains. We consider only the short spin chains because the problem of strong entanglement 
between all nodes is unresolved even for small quantum systems. 

Author thanks Professor E.B.Fel'dman for useful discussions. This work is supported by 
the Program of the Presidium of RAS No. 7 "Development of methods of obtaining chemical 
compounds and creation of new materials" . 
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